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Abstract 

We introduce a concept for random tilings which, comprising the conventional 
one, is also applicable to tiling ensembles without height representation. In particu- 
lar, we focus on the random tiling entropy as a function of the tile densities. In this 
context, and under rather mild assumptions, we prove a generalization of the first 
random tiling hypothesis which connects the maximum of the entropy with the sym- 
metry of the ensemble. Explicit examples are obtained through the re-interpretation 
of several exactly solvable models. This also leads to a counterexample to the ana- 
logue of the second random tiling hypothesis about the form of the entropy function 
near its maximum. 



1 Introduction 



Perfect quasiperiodic tilings provide useful tools for the description of certain thermody- 
namically stable quasicrystals [^, |32|, f|9| , at least for their averaged structure. However, 
it is widely accepted that the perfect structure is an idealization, and that a systematic 
treatment of defectiveness is necessary for all but a few quasicrystals |12, 18|, |30fl . In 
particular, starting from the perfect structure, the inclusion of local defects fl2| is helpful 
to overcome a number of pertinent problems, ranging from the observed degree of im- 
perfections Jl5| to the impossibility of strictly local growth rules for perfect quasiperiodic 
patterns [ pT|j . 

An alternative approach relies on a stochastic picture and adds an entropic side to the 
stabilization mechanism. Its most prominent representatives are usually formulated in 



terms of so-called random tilings [ITR E9|, |30fl . It has long been argued that such stochastic 
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tilings should provide a more realistic description both of the structures with their intrinsic 
imperfections and of their stabilization, which is then mainly entropic in nature. 

Recently, improved methods have been suggested fl34| to test this hypothesis in prac- 



tice, and the results support the claim that such models are realistic and relevant [35]. 
Further impetus is given by the very recent progress to grow (almost) equilibrium random 
tilings on the computer J33] in suitable grand-canonical ensembles. 

The concept of entropic stabilization applies to a wide class of - not necessarily qua- 
sicrystalline - tiling ensembles (defined through a set of proto-tiles and packing rules). 
However, within the original approach to random tilings (see |30 for a review), only a 



(while interesting) rather special subset of these can be described. This is because two 
concepts are mixed here that, in principle, have nothing to do with each other: the mech- 
anism of entropic stabilization and the concept of the phason strain, the latter being used 
for the description of quasicrystalline order for those tilings that allow a so-called height 
representation. Hereby, the phasonic language is used throughout, restricting the range 
of applicability. 

In this paper, we propose a more elementary and, at the same time, more rigorous point 
of view which directly employs the densities of the tiles to describe the phase diagram. 
While allowing for adequate generality, the new picture conforms to the other one, if this 
exists - at least, if the relation between the tile densities and the strain parameters is 
locally linear. 

This reformulation is, however, not done for mere aesthetic reasons, but has a number 
of advantages and leads to consequences of direct interest. In particular, the generality 
of the approach now makes classes of exactly solvable models accessible - e. g. examples 
compatible with crystallographic symmetries - which already show a variety of interesting 
and typical phenomena. 

For an explicit result, let us mention the two so-called random tiling hypotheses. The 
validity of these hypotheses is the usual set-up for an elastic theory in order to show the 
existence of a diffraction pattern with Bragg peaks for three-dimensional tilings, see 



and, for a recent discussion of experimental results, |I7| . The first hypothesis assumes 
that the point of maximum entropy in the phase space is one of maximal symmetry. We 
will show below how it can be derived from more basic assumptions as a consequence, 
and need not be stated as hypothesis. The second hypothesis uses the phasonic picture of 
quasicrystals to establish a kind of elastic theory. The hypothesis would now state that 
the entropy is locally a quadratic function of the densities near its maximal value, and the 
Hessian can be interpreted as an entropic elasticity tensor. However, this picture should 
to be taken with a pinch of salt: we shall show a counterexample where the maximum is 
unique, but coincides with a second order phase transition and is thus not of quadratic 
nature. 

In the conventional setup much information, such as elastic constants and diffraction 
properties, are extracted from the correlations of the height variables. The fundamental 
variable analogous to the height correlations in our picture turns out to be the covariance 
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matrix of the tile numbers. 

Let us sketch how the paper is organized. The following section is devoted to the 
introduction of our concept and the thermodynamic formalism used. After a brief discus- 
sion of the class of (generalized) polyomino tilings p3| , which comprises all our examples, 



we will illustrate our concepts by simple examples in one dimension. In the section on 
symmetry versus entropy, we present an argument on how to derive the first random tiling 
hypothesis from elementary assumptions. This is then substantiated by various planar 
examples, where our focus is on the re-interpretation of exactly solvable models which 
thereby gain a slightly different interpretation and, eventually, even another application. 
Our examples are taken from dimer models, but also the three-colouring model on the 
square lattice and hard hexagons are discussed. 

2 Setup and preliminaries 

In this section, we describe the way how we deal with random tilings. Since we do not 
want to use the special and somewhat restrictive phasonic picture for our fundamental 
definitions, these will differ from those given elsewhere [|30|j . 

As a (random) tiling we define a face-to-face space filling with tiles from a finite set 
of prototiles, without any gaps or overlaps. There might be a number of additional 
local packing rules specifying the allowed tilings. The number of allowed patches should 
increase exponentially with the volume of the patch (and hence with the number of tiles), 
resulting in a positive entropy density for the tiling ensemble, which we then call random 
tiling ensemble. For the tiling ensemble being homogeneous, prototile shapes and packing 
rules should admit tilings containing all sorts of tiles (with positive density), otherwise the 
ensemble splits into different subensembles that have to be taken into account separately. 



With this setting, random tilings in the sense of |3(| are included as a proper subset. 
In order to give the thermodynamic limit a precise meaning, we formalize the definition 
as follows. 

For a random tiling of the space R d , its prototiles % are compact, connected subsets 
of R d of positive volume k, homeomorphic to balls. Let A C R d be given as connected, 
compact set of positive volume V(A). We now have to explain our notion of a patch. If 
Pa is a collection of translated prototiles, we call Pa a A-patch iff Pa is a covering of A 
(without gaps) such that each element of Pa has a nonempty intersection with A, and 
every two elements have empty common interior. We call two A-patches equivalent iff 
they are translates of each other. For fixed prototile numbers n>i, . . . ,um (with corre- 
sponding densities = lirii/V(A)), let us denote the number of non-equivalent A-patches 
by gtSjhi, . . . ,n M ). 

It should be stressed here that, in principle, we can impose any type of boundary 
conditions for our definitions. However, free boundary conditions are the natural ones to 
choose from a physical point of view - a different choice may indeed lead to deviating 
results; fixed boundary conditions in particular seem to be too restrictive in most cases 
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|]54l , |Tg| ? [43f . On the other hand, most of the exact results have been obtained by use of 
periodic boundary conditions, and it has to be shown, for each example, that the results 
also apply to the case of free boundary conditions. 

We will use a grand-canonical formulation in the following, in difference to the (usual) 
canonical treatment of tiling ensembles [[$(], 4"1| , as we believe that an ensemble with fluc- 



tuating tile numbers is the adequate setup to cover processes such as local rearrangements 
or tiling growth. Moreover, this formulation will prove advantageous for our arguments 
and practical calculations. 

As we are interested in configurational entropy only, we will assign equal (zero) energy 
to every prototile. In this situation, the chemical potentials of the prototiles are the only 
parameters remaining. As variables conjugated to the densities, we introduce chemical 
potentials /ij as well as activities Zj, = e Mi for every prototile, Q (i = 1, . . . , M), and define 
the grand-canonical partition function to be the configuration generating function 



Z A (//i,..., //m) = E 9A(n 1 ,...,n M )z% rh ---z$ nM . (1) 

ni,— ,um 



Since we are interested in the entropic behaviour of large systems, we define the grand- 
canonical potential Q per unit volume to be the limit 

0(/ii,...,/i M ) = lim — —log Z A (//i,..., Vm)- (2) 

A^oo y \A) 

Here, {A} is any collection of connected, compact subsets of R d of positive volume V(A) 
such that lim\^(A) = oo and lim5(A)/V r (A) = 0, where -B(A) denotes (the maximum) 
over the volume of the boundary tiles of any A-patch, in the same spirit as the van Hove 
limit taken for ensembles in statistical mechanics ||55|| . For models of physical relevance 
the thermodynamic limit defined that way should exist for all physical quantities and be 
independent of the special choice of the limit sequence {A}. 

The mean (volume) densities of the different prototiles in this ensemble can be com- 
puted as 

Pi{p-x, ■••,Hm)= . lim 77TTs(kni)A = ^ ^' ' " ,m (i = 1,...,M), (3) 

A^oo V (A) Ofli 

where ( )a denotes the finite-size ensemble average for given chemical potentials //i, . . . , [Im- 
Note that all physical quantities defined here are ensemble averages. In cases where 
these quantities differ from corresponding values of a typical representative of the ensem- 
ble, one has to reflect whether the grand-canonical setup fits the experimental situation, 
or if a canonical ensemble is the correct choice. Bear in mind, however, that the van Hove 



1 We suppress extra prefactors /? in all definitions below. 

2 If Boltzmann-factors are taken instead of activities, this corresponds to the (canonical) free energy 
in statistical mechanics. 
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limit is not properly defined for single tilings that are not self-averaging. Thus either 
canonical and grand-canonical thermodynamics coincide, or the the canonical picture is 
likely not to be well defined at all. 

So far, the chemical potentials and also the conjugated (mean) densities do not form 
an independent set of macroscopical parameters to describe the ensemble. We always have 
T,iL\ Pi = 1 as normalization constraint, and in general, the tile geometries and packing 
rules may eventually impose further constraints on the densities. 

In order to make the notion of macroscopical independence precise, let v = v j n j 7^ 
be an arbitrary linear combination of some prototile numbers nj, with real coefficients v j. 
We call the accompanying set of tile densities {pj} independent in the thermodynamic 
limit, iff the variance per unit volume of v is strictly positive for any choice of the Vf. 

^<(«-<*>a) 2 >a > 0. (4) 

The degrees of freedom of the tiling ensemble are given by a maximal set of independent 
densities. One checks that, with this definition, normalization reduces the independent 
densities by one, and there is no additional freedom due to surface effects or subclasses 
of the ensemble that do not contribute entropically. 

In general, since the average taken in (ffl) depends on the chemical potentials, the 
variance of some v may vanish only for a special choice of the potentials leading to a 
reduction of independent densities at isolated points of the phase space. However, no 
such thing will happen in any of the examples we discuss in this sequel and, for simplicity 
of the presentation, we exclude this possibility in the following. 

Let us mention a consequence of this definition for the covariance matrix of the tile 
numbers (which corresponds to the isothermal compressibility in the thermodynamic sit- 
uation), in the thermodynamic limit: 

Kij = li^y^^hniljU^A - {kni)A{ljnj)xJ = (i,j = l,..., M). (5) 

This matrix describes tile fluctuations around the mean density and correlations among 
different tiles. A straightforward calculation shows that ([5]) is indeed strictly positive iff 
it is restricted to the subspace of density parameters that had been chosen independently 
according to the above definition. 

In general, a linear constraint on the tile numbers of the form 

M 

VjiljUj) = cV(A) + 0(V a (A)) (a < 1, c 6 R) (6) 

leads to linear variation of the grand-canonical potential in certain coordinate directions: 
0(^i+ Wit, ...,n M + v M t) = (f>(fii, . . . , fj, M ) + ct (teR). (7) 
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For each constraint of this form, we can shift a single chemical potential to zero in the 
grand-canonical potential by a suitable choice of t. 

Let us choose an independent subset of macroscopic parameters, p~\, . . . , pk say, which 
we abbreviate as {/?}, and write the random tiling entropy as a function of these indepen- 
dent densities. According to the considerations above, this can be done as follows. We 
set fik+i — • • ■ = A*m — i n (E)) i- e - we normalize the grand-canonical potential by an ap- 
propriate choice of the origin of the energy scale. We call the normalized grand-canonical 
potential <ft again. From (0) and (^) it is readily verified that the normalization of <fi leaves 
the independent densities invariant. 

The (grand-canonical) entropy per unit volume as a function of the independent 
(mean) tile densities is obtained as Legendre transform of the (normalized) grand-canonical 
potential [] with respect to the conjugated chemical potentials, 

k 

s(pi,...,p k ) = <f>(pi({p}),...,Hk({p})) - ^PiPi{{P})- (8) 

i=i 

Note that the Legendre transform is well defined owing to the non-singularity of the 
matrix 

dp % dp 3 ~ KtJ iM-i,...,*;- w 

Let us locate the maximum of the entropy curve as a function of the independent 
densities pi, . . . ,p k . We have by definition 

= -&{{?}) (; = !,...,*). (io) 

It follows that the entropy is extremal, s max = 0, at p% = ■ • • = pk — 0. (For the grand- 
canonical potential @ this point occurs at p± = ■ ■ ■ = pm — 0.) This is the situation of 
the maximal random ensemble where each configuration is weighted equally. Let us take 
a closer look at the Hessian, the matrix of the second derivative 



dh({p}) _ dp 3 {{p}) 



«tm. (ij = !,...,*). (n; 



dpidpj dpi v J v 

If s is locally C 3 , we conclude that the entropy is a concave function everywhere. On 
the other hand, also in the vector space of the independent densities, the variance of v 
and hence the covariance matrix (|5|) may diverge in the thermodynamic limit for certain 
values of the chemical potentials. This generally occurs at phase transitions. In this case 
the Hessian of the entropy has as an eigenvalue in the large system limit; we will give 
an example below where this happens at the point of maximum entropy. For all other 

3 Under the above assumption, this definition coincides with the micro-canonical one ; for a number 
of examples in the sequel this can be shown explicitly. 
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situations, the shape of the entropy curve near its maximum is completely determined by 
the Hessian. 

We interpret the negative of the Hessian, at the point of maximum entropy, as tensor of 
entropic elasticity. Without any further symmetry, we refer to its eigenvalues as the elastic 
constants of the random tiling ensemble. We will define a symmetry-adapted version in 
the section on symmetry. 



3 Polyomino tiling ensembles 

A class of (crystallographic) tilings that allow a re-interpretation of many exactly solved 
models are the (generalized) polyomino tilings, see [23, [|3j for definitions and background 
material. Their prototiles are taken from one or combinations of several elementary cells of 
a given periodic graph. There may be as well certain packing rules, resulting in restricted 
random tiling ensembles. We will give a number of examples of tilings consisting of 
monominoes and dominoes. 

The simplest class of random tilings are tilings of (coloured) monominoes without 



further restrictions. For lattices of any dimension, this results in Bernoulli ensembles |53 
A natural way to impose packing rules for these models is to exclude the occurrence of 
clusters with moniminoes of the same colour. For nearest-neighbour exclusion on lattices, 
the restriction is hierarchical: the tiling can be regarded as being built from consecutive 
layers of "hypertilings" with additional constraints along the layer direction. In this case, 
the entropy is bounded by the entropy of the corresponding hypertiling ensemble. Several 
two-dimensional examples are discussed below. 

We should mention the one-to-one correspondence between polyomino tilings on a 
periodic graph A and polymer arrangements on the dual cell complex A* (the so-called 
Delone complex). This is illustrated in figure [l] for a simple case, the monomer-dimer 
system on the square lattice. 



Figure 1: Each polyomino tiling is a polymer configuration on the dual cell complex. 



Owing to this connection, the (pure) domino problem can be solved in a number of 
cases in two dimensions, as the dual fully packed dimer problem allows the exact solution 
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on planar graphs by means of Pfaffians Q. This has been shown by Kasteleyn in [[38], [3£ 
where he computed the entropy of dimers on the square lattice. It is even possible to 
compute densities of finite patches explicitely, which has recently been shown by Kenyon 
40fl . We will treat two examples representing the different types of critical behaviour of 



the solvable domino systems ||48|| : systems of Onsager type (with logarithmic singularity 
of the covariance matrix) and systems of Kasteleyn type (with square-root singularity of 
the covariance matrix). 

Let us consider monomino-domino tilings without restriction next. The one-dimensional 
system will be discussed below. In two dimensions, there are various connections with 



other models, e. g. for the triangle-rhombus tiling on the triangular lattice [p9| , but, 
to our knowledge, no explicit exact solutions. On the other hand, as follows from the 
theory of monomer-dimer systems p8|| , phase transitions can only occur for tilings with 
vanishing monomino density, regardless of the dimensionality of the system - these phases 
are the (pure) domino tilings mentioned above. There is at least one monomino-domino 
tiling with constraints that allows for exact solution in two dimensions: the tiling on the 
hexagonal lattice together with exclusion of neighbouring dominoes on different sublat- 



tices. This tiling ensemble is equivalent to the ensemble of hard hexagons |9], [11]], see also 
below. 



4 Unrestricted and one-dimensional examples 

The simplest class of random tilings, built from M coloured monominoes on an arbitrary 
lattice, can be described purely combinatorially. There are no local constraints or packing 
rules, and thus no global constraints for the densities besides normalization. One applica- 
tion of this model is the description of unrestricted chemical disorder in multicomponent 
alloys §. 

The thermodynamic limit may be taken by counting the number of monomino config- 
urations on consecutive lines of length N. In this simple case, the densities pi [] can be 
given as explicit expressions of the activities: 

P*= 7 Zj — (i = l,...,M). (12) 

Z\ + . . . + Z M 

In order to eliminate the normalization constraint, we choose p±, . . . , pu-i as independent 
densities and normalize the grand-canonical potential via 

0(/ii, . . . , Pm-i) ■= <f>(fJ>i, Pm-i, 0). (13) 
Both functions are related via ([?[), 

4> {Hi - Hm, Hm-i - Pm) = 0(^1, • • • , Pm) - Pm- (14) 

4 For an introduction into the Pfaffian method, we refer to Appendix E of |57j]. 

5 We suppress the overbar and write pi instead of pi from now on. 



S 



We obtain the entropy as Legendre transform of the normalized grand-canonical potential, 



M 

s(p u . . . , Pm-i) = -^PjlogPi, (15) 

3=1 

where pu = 1 — Hf=\ l Pj- This coincides with the entropy of a Bernoulli system |xj with 
fixed frequencies pi = p%. The entropy function is strictly concave and takes its (unique) 
maximum at the point pi — ■ ■ ■ — Pm = -a, with value 

M 1 1 

W = -E M log M =logM - ( 16 ) 

3=1 

Near its maximum, the entropy curve is a quadratic function of the independent order 
parameters. Its shape near this point is completely determined by the Hessian 



d 2 s 



(Pm+PT%) . =-M(l + 5 ZJ ) (*,j = l,...,M-l). (17) 



dpidpj 

The negative Hessian is the elastic tensor. In particular, for M = 2, we obtain 



S (p) = log2-2(p- 2 




which (with E — 2p — 1) conforms to the phasonic expression given in [|30[ 



Owing to the lack of geometric or other constraints, this class of examples is fully 
controlled by the strong law of large numbers 0: With probability one, for a given 
member of the ensemble, the probability to occupy a position with a tile of type j is its 
frequency, pj. This even connects the growth with the equilibrium properties, as also 
observed recently for other, more complicated ensembles [33| . 



A first natural step to examine the effect of different tile geometries as well as of 
certain packing rules is to allow these in one coordinate direction. Introducing different 
tile sizes, one of the simplest models is the ID monomino- domino model PS . 



Similar to the previous example, the only constraint of this model is due to density 
normalization. We eliminate this constraint in the same way, i.e. we write 

N 

Z N (p) = Y,9 N (k)z k (19) 

fc=0 

where gN{k) is the number of connected patches with iV occupied sites, k of which being 
monominoes. So, p = log z is the chemical potential of a monomino. Obviously, Z^{p) 
fulfils the Fibonacci-type recursion formula 

Z N (p) = zZ N _ x {p) + Z N _ 2 (p) (20) 
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Zq : = 1 Z\{n) = z. 

The quotient Z^{p) j Z^-\{p) is a rational function whose limit, for iV — > oo, exists and 
can easily be calculated from (|20|) to be 



lim^U IZ + Ji+fE)*). (21) 

The grand-canonical potential per site is defined through <j)(p) = liniAr^oo ■= \ogZ N (/i), 
but can directly be calculated from the quotient (]2lf) as 



0(p) = lim log = arsinh^. (22) 

We denote the monomino density by pi, the domino density by p 2 (pi + p 2 = 1). Let 
us now write the entropy curve as a function of the monomino density p±. The entropy 
curve can be computed from the grand-canonical potential as 

s(pi) = 4>{z(px)) -pilog^(pi) (23) 



(Pi + y)l°g(Pi + y) - Pi log Pi - y lo Sy' 



Pi W = — = tanh0 p 

<9p / x / . ^ 



where 



> 2 y 

The entropy function is a strictly concave function with a quadratic maximum 

s m ax = logr^ 0.4812 (24) 

a ^ ^pt = T ~ ^-6180) where r = 1+ 2 V ^ denotes the golden number []. We obtain an elastic 
constant A = \Vh « 2.7951. 

It is possible to derive a closed expression for the grand-canonical partition sum |28[ . 
Defining Pn{x) = i~ N Zn(2ix) transforms the recursion (|20|) into that of the Chebyshev 
polynomials of the second kind |l[ which are 

P N(cos()) ^ME±ni (25) 
sin 

wherefrom one gets 



even 
odd 



(26) 



6 The silver number A = 1 + \/2 appears in the monomino-domino model of monominoes with two 
different colours. 
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Packing rules and supertiles 

So far, no further constraint or packing rule has been imposed on the ensembles. A 
natural way to do this in one dimension is to exclude certain strings of consecuting tiles. 
In many cases, these ensembles can be directly dealt with by recursive methods as shown 
above (e.g. the M-colouring problem in ID), or the constraints can be eliminated by 
a reformulation of the problem in terms of 'supertiles'. Let us illustrate this for the 
ensemble of chains in two monominoes a, b with the exclusion of r consecuting a's and s 
consecuting 6's. A moments reflection shows that any allowed chain can (up to boundary 
tiles) uniquely be written in supertiles A it j, 1 < i < r, 1 < j < s, forming strings of i 
consecuting a's, followed by j o's. Vice versa, any sequence in the supertiles translates to 
an allowed one in a,b. 

Similarly, the ensemble of monomino chains with tiles a, b and exclusion of bb is equiv- 
alent to the monomino-domino model mentioned above. We give the entropy for this 
ensemble via its 'supertile' formulation with polyominoes A = a and B = ab, starting 
from a Bernoulli ensemble rather than taking different tile geometries of the supertiles 
into account. The frequencies of the two models are related according to 

Pa+Pb Pb 

Pa = -z — , Pb = — : — , 27 

p A + 2p B p A + 2p B 

, Pb Pb 
Pa= 1 , Pb = —■ 

Pa Pa 

It is obvious from here that the possible values of p a ,Pb are restricted to 

\ < Pa, 0< Pb < Pa <l. (28) 
As the entropy per supertile is given by 

s(pa,Pb) = -Pa log Pa - Pb \ogp B , (29) 
the entropy per (small) monomino follows as 

s(Pa,Pb) = (l + — I s(p A (p a ,Pb),PB(Pa,Pb)) 
V Pa) 

= p a logp a - p b logp b - {p a ~ Pb ) \og(p a - Pb ). (30) 

Taking x = 2p a — 1 as parameter, < x < 1, the entropy curve is exactly the curve of 
the monomino-domino model. The maximum occurs at p a = ~ 0.7326. We obtain 
an elastic constant A = 5V5 11.1803. 

However, already in ID, more complex examples may be considered. An interesting 
family is provided by the ensemble of square-free words in n letters, n > 3, which is known 
to display positive entropy. In fact, the value 

s(n) = arcosh ( - (31) 
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gives (for n > 4) a reasonable lower bound of it, which is asymptotically exact, but the 
determination of the exact value is still an open problem 0. 

5 Entropy and symmetry 

In this section we consider effects of the symmetries of a tiling ensemble on its entropy 
function. It should be stressed again that the results below are derived using the grand- 
canonical formalism. All physical quantities are averages over the whole tiling ensemble. 
In most situations of physical relevance, however, the ensemble averages coincide with the 
corresponding quantities of a typical single tiling, and all results translate to the canonical 
ensemble as well. This will be the case for all examples discussed in this article. 

As a symmetry of a random tiling ensemble, we define every linear bijection S on the 
density parameters that leaves the entropy function invariant. We write 

M 

•V/V •/;, E<Vi (i = l,...,M). (32) 
3=1 

Owing to the normalization constraint, the restriction of S p to the subspace of independent 
density parameters where the entropy function is defined leads in general to an affine 
bijection. 

Examples are obtained from bijections on the tiling ensemble which transform den- 



sities like (32). They are symmetries of the random tiling ensemble if they respect the 
thermodynamic limit, i.e. if they induce a bijection on the set of A-patches, for each A 
of a suitably chosen limit sequence {A}. The simplest examples are colour symmetries, 
which permute equal prototiles of different colour. Moreover, rotations and reflections 
on the tiling ensemble, together with a limit sequence of circular sets, lead to geometric 
symmetries. Let us denote every symmetry which is not of the former type as hidden 
symmetry. We will discuss a random tiling with a hidden symmetry below. 

In general, every linear bijection on the densities can as well be represented in the 
vector space of all chemical potentials; as can be verified from ([3]), a transformation S p 
induces the adjoint mapping 



M 



S„ : /i« -> fit = E ($% N (< = 1, . . . , M), (33) 

3=1 

where S* = {S^ 1 ) 1 means the matrix inverse and transpose of S. In particular, the imparl- 
ance of the grand-canonical potential under symmetry operations ([53]) imposes symmetries 
of the (mean) density vector field in the space of chemical potentials, 

p(fi)=p((i). (34) 

We conclude immediately that symmetry related tiles having equal chemical potential 
occur with the same (mean) density. Note that the point of maximum entropy (fi = 0) is 
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a fixed point under any given symmetry which means that symmetries constrain possible 
density parameters at the entropy maximum. This is precisely the statement of the so- 
called first random tiling hypothesis: 

• The point of maximum entropy is a point of maximum symmetry. 

For random tilings with a height representation, maximum symmetry implies zero phason 
strain, E = |30| , so this is contained. In our setting, the conlusion above implies: 



At the point of maximum entropy, symmetry related tiles occur with equal (mean) 
density. 



We can further conclude, as can be read off from (p4|), the density vector p, at the point 
of maximum entropy, has only components in direction of the trivial one- dimensional 
representations of the symmetry group. Thus, the symmetry determines the entropy 
maximum completely, if the representation of the symmetry group on the vector space of 
independent density parameters is irreducible (and nontrivial). We will give examples of 
that kind later on. 

Let us now focus on the transformation behaviour of the entropy function (RI). This 
is done most easily within a symmetry-adapted parametrization. By an affine coordinate 
transformation, it is always possible to introduce new independent density parameters 
such that the new origin is the point of maximum entropy, and the symmetries are repre- 
sented by linear, orthogonal transformations. In many cases, the new density parameters 
have a geometric interpretation in form of densities of supertiles. Let us call these new 
parameters r and assume that we are in such a frame. 

The entropy is invariant under symmetry transformations, 

s(r) = s(f). (35) 

At its maximum, the entropy expansion can be written in terms of invariants of the 
given symmetry group |_]. The form of the second order terms is determined by the 
Hessian H of the entropy (the matrix of its second derivatives). The symmetries induce 
a transformation behaviour of the form 

H\ r = (S r Y H\ f S r . (36) 

Since S r is an orthogonal transformation, it commutes with H at the point of maximum 
entropy, 

H\ S r = S r H\ . (37) 

According to Schur's lemma [27], H\ Q acts trivially on the irreducible subspaces of the 
symmetry group and can thus be written as a linear combination of projectors P^ 1 onto 
the irreducible components, 

#| = -£A»P« (38) 

7 For the quasicrystal entropy as a function of perp strain, this kind of arguing was used in M] . 
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We refer to the Aj as elastic constants of the given random tiling. This is more appropriate 
than the use of the eigenvalues since they are no longer independent in the presence of 
symmetries. The quadratic term of the entropy expansion is then of the form 

s 2 (r,r) = -i5>/«(r,r), (39) 

where the jW(r, r) are the quadratic invariants defined by 

/W(r,r)= <r,P w r>, (40) 

and < •, • > denotes the scalar product. We will give examples later on. 

Now, the discussion on the relation of entropy and symmetry can be extended even 
further. Above, we have always taken a set of tiles as basic elements to construct the 
tiling. This is, however, not the only possibility. Instead, tilings can also be seen as 
built from bonds, vertex configurations, patches of a given volume or number of tiles, or 
the like. Adopting this point of view, we can assign chemical potentials and densities to 
patches or bonds rather than to the prototiles. While the number of basic elements and 
packing rules may increase quite rapidly this way, the above arguments can nevertheless 
be repeated without change, and we conclude that ensemble densities of symmetry related 
patches of any size are equal at the point of maximum entropy. If the density parameters 
of a typical tiling of the ensemble coincide with their ensemble averages, then: 

• Entropically stabilized random tilings are (with probability one) locally indistin- 
guishable^ from their images under symmetry transformations. 

This is the main result of this section. Note that, also for tilings with a height represen- 
tation, this is a much stronger statement than the demand E = on the phason strain, 
since different densities of symmetry related patches may be compatible with the same 
value for the phason strain. 



6 Examples in two dimensions 



6.1 



The rhombus tiling 
representation 



a domino tiling of Kasteleyn type with height 



The rhombus random tiling [19], |Tj] is the ensemble of all coverings of the plane with 60°- 
rhombi. It is also called a lozenge tiling. It is a domino tiling on the triangular lattice, 
equivalent to the fully packed dimer model on the hexagonal lattice [3S, 5S]. We discuss it 
here since it is the prototype of a random tiling possessing a height representation, which 
allows us to demonstrate the equivalence of our approach to the conventional one, for 
this model. In the following, we focus on the symmetries of the tiling and give a simple 



8 An introduction to equivalence concepts of tilings is given in 0] 
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derivation of the entropy. In the appendix, we relate the elastic constant given by Henley 



29], |30j to our result. 

The symmetry of the random tiling ensemble which transforms rhombus tilings into 
other allowed ones is Dq = D% x C 2 . For undecorated tiles, however, the inversion sym- 
metry of the prototiles already enforces C 2 as (minimal) symmetry for any element of 
the tiling ensemble, thus 5*3 ~ D 3 remains as relevant geometric symmetry in the space 
of density parameters. There are no hidden symmetries in this case. Since all pro- 
totiles are related by rotation, for the point of maximal entropy we immediately conclude 
Pi = P2 = P3 = §• 

The filling constraint pi + p 2 + p 3 = 1 reduces the number of independent variables to 
two. As the symmetry acts irreducibly on the remaining two-dimensional subspace, the 
symmetry of the entropy function must belong to a two-dimensional irreducible represen- 
tation of -D3. It is most convenient to represent the phase space as the set of points inside 
an equilateral triangle, see figure |2], this way obtaining a handy parametrization which 
we call its phase diagram. We choose its center as the origin and the vectors pointing 




r l r 2 

Figure 2: The phase space of the rhombus tiling is an equilateral triangle. 



to the corners of the triangle as unit vectors of the different tile densities. The center 
corresponds to the point of maximal entropy, and each point r = (ri,r 2 ,r 3 ) corresponds 
to an ensemble with densities pi = r« + |. The symmetries of the entropy function are 
then reflected in the symmetries of the triangle. The point of maximum entropy is the 
only one with D% symmetry. Along lines with equal density of two different tiles there is 
reflection symmetry. 

The grand-canonical potential per rhombus as a function of the different tile potentials 
Pi = logZj (with i = 1, 2, 3) reads |58[ 

2vr 2tt 

0(/ii,/i 2 ,/i 3 ) = J J logdet X((p u (p 2 )d(p 1 dcp 2 , (41) 




2- 



— J log (max , z\ + z\ - 2z 2 z 3 cos dip, 
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where the determinant is given by 

det A(<pi, <p 2 ) = z\ + z\ + z\ + 2zi z 2 cosy?i + 2^2:3 cos <p 2 + 2zsZi COs(<Pi - <p 2 ). 

This result was obtained by taking rhombohedral patches of increasing size, together 
with periodic boundary conditions. The result applies nevertheless to the corresponding 
limit with free boundary conditions, as can be seen as follows. The rhombus tiling can 
be regarded as a limiting case of the dart-rhombus tiling described below, where some 
chemical potentials are set to — oo. This operation commutes with the free thermodynamic 
limit [46], yielding the same result as in fl4"I]) . 

We introduce the normalization constraint by setting 

0(/ii,/i 2 ) := 0(/ii,/! 2 ,O). (42) 
The relation between <j) and <j) is given by 

4> (Pi _ ^3, P2 - A*s) = <f>(Pi, P2, P3) - A*3, (43) 

Pi (Pi - A*3, ^2 - ^3) = Pl(/il, /i2, P3), P2 (Pi - P3, P2 ~ Pz) = Pz(Pl, P2, P3) ■ 

Let us determine the entropy curve as a function of the independent densities p\ = p\ 
and p 2 = p 2 . From the normalized grand-canonical potential <ft we obtain 

f- ~\ sin 71 Pi ( ~ ~n sin np2 

Zi{pi,p2) = - , _ , rr , Wi,P2) = - ,~ , - v 44 

sin 7r(pi + p 2 ) sin 7r(pi + p 2 ) 

The expression for the entropy follows as 

s(pi,P2) = - log z 1 (p,p 2 )dp = - / log- ———dp. (45) 

io jo sm7r(p + p 2 ) 

At the point of maximum entropy we find 



1 ~ sin^fc 

W = - E ~rf- ~ °- 3230 > ( 46 ) 
n k=i fc 

which is considerably smaller than the value ^ J2kLo (2fc+i) 2 ~ 0.583 of the square lattice 
domino problem [ j3Sfl . 

As the irreducible representation of the symmetry group is two-dimensional, we expect 
to have a single elastic constant. In fact, for the quadratic term in the entropy expansion 
we find 

1 71 

s 2 {r,r) = --■ -j= + m + (47) 
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On the other hand, the entropy function is defined on the subspace complementary to the 
vector (1, 1, 1)*. This yields 



J(r, r) = 2(r\ + r 1 r 2 + r|). 



(4f 



Thus we find an elastic constant of A = ^j=. In the appendix we relate this constant to 
the constant given by Henley [29, 30 1 which was defined from the height representation 
of this model. 

Note that there is a phase transition at points where two densities approach zero. In 
the above parametrization, this corresponds to the corners of the triangle. At these points, 
the covariance matrix shows a square- root singularity [|14j], indicating a phase transition 



of Kasteleyn type |48[ . 

We give an expression of the entropy along a line of reflection symmetry: in the regime 
Pi = P2, we have 



7T JO 



7rpi 



log 2 cos xdx. 



This function is plotted in figure |3|. The elastic constant is A = 



(49) 



0.3 
0.25 

0.2 
0.15 

0.1 
0.05 





Figure 3: Entropy of the rhombus random tiling in the regime px = Pi- 



6.2 The dart-rhombus tiling — a domino tiling of Onsager type with hidden 
symmetry 

We now present another random tiling of the domino class, but this time with two different 
types of prototiles Q. It is therefore slightly more interesting than the model discussed 
above. The prototiles of the dart-rhombus random tiling are 60°-rhombi and darts made 
from two rhombus halfs, see figure |[ In addition to the usual face-to-face condition, 
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we impose an alternation condition on the rhombi, such that neighbouring rhombi of 
equal orientation are excluded. Finally, to avoid pathologic lines of alternating darts, we 
demand that two neighbouring darts must not share a short edge. These tiling rules force 



the darts to arrange in closed lines, see figures |10| and |TT[ We thus deal with a dilute loop 
model of dart loops in a background of rhombi. For more on this picture, which we will 
not expand here, see [j56 ]. The dense loop phase has the minimal total rhombus density 
p = \, see also figure |[ 

This random tiling corresponds to the fully packed dimer system on one of the so- 
called Archimedian tilings (where all faces are regular polygons) . Lattices and tilings 



of this type have already been described in detail by Kepler |^TJ, but have since been 
rediscovered many times. In the context of statistical mechanics the graph of our model 
has also been called a Fisher lattice []57j, its dual graph is named a Kagome lattice. The 





Figure 4: Prototiles of the dart-rhombus tiling and fundamental cell of the 
Fisher lattice. 



first solution of the Fisher lattice dimer model was given by Fan and Wu |2(J by realizing 
the correspondence to a soluble subcase of the the eight- vertex model fTI], the so-called 



free-fermion case. A one-to-one mapping between the free-fermion case of the eight- vertex 
model | and the dart-rhombus tiling is obtained by a reformulation of the latter in terms 
of 'supertiles', see figure ||. Closed dart loops in the tiling correspond to closed edge 
loops on the square lattice. 

Another interpretation exists as a model with domain walls where different walls can 
meet and annihilate 



13 



in generalization of the dimer system on the hexagonal lattice, 
which can also be interpreted as a domain-wall model. In our model, the domain walls 
correspond to the dart loops. 

Let us denote the densities of rhombi by pi (i — 1, 2, 3) and the densities of darts by 
Oi (i — 1, ... ,6). We have several constraints on these macroscopic variables: in addition 
to the filling condition, the loop condition means that densities of opposite darts have to 



9 We follow the Baxter's convention here pj] ]. 
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Wi = yiz^ w 2 =y!y 2 5§ w 3 =y 3 z^ w^y^z. 




Figure 5: A mapping between the eight-vertex model and the dart-rhombus tiling via "super- 
tiles". 



be the same, 

o"i = cr 4 , a 2 = 05, 0-3 = cr 6 . (50) 

Moreover, as each dart is accompanied by a rhombus of equal type, the remaining rhombi 
occur with equal frequency of each type due to the alternation condition, 

Pl - ^ = p 2 - cr 2 = p 3 - 0-3. (51) 

These constraints reduce the number of independent parameters to three. 

The group of geometric symmetries of the tiling is D%. But these are not the only 
symmetries of the model: there is a hidden symmetry which exchanges darts and rhombi. 
It is most easily described in the vertex model formulation where it corresponds to the 
exchange of free and occupied edges (spin reversal symmetry). 

The full symmetry group of the dart-rhombus tiling turns out to be the tetrahedral 
group Td- As the group action respects the density constraints and acts irreducible on 
the remaining subspace of independent density parameters, we can characterize the phase 
space of the dart-rhombus random tiling according to the three-dimensional irreducible 
representation of the tetrahedral group T^. Each phase of the dart-rhombus tiling cor- 
responds to a point inside the tetrahedron: we fix the tetrahedron center as the origin 
and choose an orthonormal basis of unit vectors pointing to edges along twofold axes, as 
is indicated in figure |6]. In this way we obtain a parametrization in terms of reduced 
rhombi densities 

n = Pi-\ (2 = 1,2,3). (52) 

The vectors pointing to the four corners of the tetrahedron are unit vectors for the reduced 
dart densities S{ and for the reduced total rhombi density r, 

Sl = 2o-i-\ {1 = 1,2,3), (53) 
r = n + r 2 + r 3 . 
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Figure 6: The phase space of the dart-rhombus tiling is a tetrahedron. 



In this geometry, the dense loop phase is found on the triangle surface with minimal 
rhombus density r = — |. The other triangle surfaces which are obtained by the exchange 
symmetry are diluted phases where the darts arrange along directed lines due to exclusion 
of one dart orientation. On the other hand, the dense loop phase is an obvious decoration 
of the rhombus tiling, suggested by figure |8|. We conclude that the surface phases are of 
Kasteleyn type with phase transitions at the corners of the tetrahedron. 

The point of maximum entropy (the tetrahedron's centre) is fixed by symmetry to 
Pi = |, (Tj = ^ where the darts and the rhombi occupy half of the tiling area each. The 
value at the entropy at its maximum 

s max = ^log2^ 0.231 (54) 

follows from a geometrical consideration: in the supertile formulation, building a rhombic 
patch by starting from a prescribed configuration on the left and on the upper boundary 
of a | sector, there are in general two different possibilities to add a new tile. This results 
in an entropy of log 2 per supertile because the choice of the boundary does not affect 
the entropy value. Observing that each supertile consists of 3 ordinary tiles by area, one 
obtains (|54j). 

In general, phase transitions occur in regions where we have infinite dart lines which 



tend to minimize the enclosed area as in percolation, see figures [10] and [11]. For a more 
concrete description, we introduce rhombus activities j/i, 2/2, 2/3, dart activities Zi, . . . ,z 6 , 
and compute the grand-canonical potential per tile via Pfaffians, yielding J57| 



2vr 2tt 

24 _ 5 / / lo S det ^1,^2)^1^2, (55) 
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where the determinant is given by 



detA(y?i,y? 2 ) = y\z[z\ + y\z\z\ + y\z\z\ + y\y\y\ 
+2y 1 y 2 (z 1 z 2 z 4 z 5 - y%z 3 z 6 ) cosy?i 
+2y 2 y 3 (z 2 z 3 z 5 z 6 - y\z x z^) cos(^i - tp 2 ) 
+2y 3 y 1 (z 1 z 3 z 4 z (i - y 2 z 2 z 5 ) cosip 2 . 

This result was obtained by using periodic boundary conditions. It applies, however, 
in the case of nonzero activities, to the free case as well, since this has been shown by 
Ruelle |55| for the more general eight-vertex model. (Note that the Kasteleyn subclass 



has already been discussed above.) 

Let us now proceed with the discussion of the critical behaviour as a function of the tile 



densities. It is of Onsager type in the generic case pOfl . At points fulfilling the criticality 
condition 

W\ + u>2 + w 3 + W4 = 2 max{w%, w 2 , w 3 , W4} (56) 

the grand-canonical potential is nonanalytic, resulting in a logarithmic divergence in the 
covariance matrix. The Wi (i = 1, . . . , 8) are the Boltzmann factors of the eight- vertex 
model whose relation to the tile activities is given in figure [5|. 

The entropy as a function of the rhombi densities indeed leads to a Hessian which 
is proportional to the identity. In the four- dimensional vector space of the reduced dart 
densities and the reduced total rhombus density, the entropy is defined on the hyperplane 
complementary to the vector (1, 1, 1, 1)*. This leads to 

I(r,r) = 2(rl + rl + rj), A = 6. (57) 

Let us focus on the behaviour of the random tiling in the sixfold symmetric phase 
where each orientation of darts, respectively rhombi, occurs with the same frequency. 
The only free parameters are the total density p of the rhombi or the total density o 
of the darts, each being equally distributed over the different orientations. The entropy 
curve as a function of the rhombus density p is plotted in figure fj]. The maximum 
is already discussed above; the elastic constant in the sixfold phase is A = 6. A phase 
transition of Onsager type occurs at p = | with the usual logarithmic singularity in the 
isothermal compressibility. A phase transition of Kasteleyn type occurs at p = 1. The 
point p = ~ is the point of maximal symmetry of the previously mentioned rhombus tiling 
(its maximum scaled by a factor of 3). 

Below are four characteristic snapshots of tilings in the sixfold symmetric phase to- 
gether with their diffraction patterns, taken at different rhombi densities pi[|. The tilings 



are periodic in the vertical direction and periodic up to a cyclical shift in the horizontal 
direction. To obtain diffraction patterns, delta-scatterers were positioned in the middle 
of every prototile. Figures p|-|TT| show contour lines of the absolute value of the scatterer's 
Fourier transform, computed in arbitrary units. 
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Figure 7: Entropy of the dart-rhombus random tiling in the sixfold symmetric phase as a 
function of the rhombus density p. 



The phase transition (at rhombi density p c = |) is directly visible in the tiling. Above 
the phase transition (p > p c ) we have just some small dart islands (figure |TT|) whereas 
below the phase transition (p < p c ) they form infinite lines as in percolation (figure |TUD . 
The phase transition is also visible in the diffraction pictures. Above the phase transition 
point, we have a small separated, sixfold symmetric background which gets connected 
below the phase transition. 

At the entropy maximum (figure , the diffuse background has maximal intensity and 
is least structured whereas at minimal rhombus density (figure §) it is still connected but 
laced up. 

We mention that the height of the diffraction spots scales with the system size, as 
follows from our numerical analysis. The diffraction pattern in the thermodynamic limit 
consists of a point part (Bragg peaks) and an absolutely continous background 0]. Al- 
though this is not the generic result one expects for 2D random tilings, it is inevitable for 
this kind of crystallographic random tiling that lives on a lattice background. 



6.3 Three-colourings of the square lattice - a monomino tiling violating the 
second random tiling hypothesis 

The three- colouring model is the ensemble of tilings with three types of monominoes such 
that no two tiles of the same type are adjacent. Since the packing rule is of nearest- 
neighbour type, the restriction is hierarchical, and it is instructive to look at the one- 
dimensional model first. In one dimension, using the recursion method demonstrated 
above, the grand-canonical potential can be shown to be the logarithm of the positive 
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Figure 8: Minimal rhombi density (p = 0.33) 



real root of the equation 

x 3 - {z\Zi + z 2 z 3 + z 3 zi)x - 2ziz 2 zz = 0, (58) 

the Zi being the activities of the different monominoes. By construction, the entropy of 
this ensemble is an upper bound for the entropy of its generalization to two dimensions. 

The group S3 of colour permutations is a symmetry for both models which also respects 
the filling constraint p± +p 2 + j°3 = 1- Consequently, the maximal entropy occurs at p\ = 
P2 = P3 = |j and the thermodynamic functions are invariant under the two-dimensional 
irreducible action of the symmetry group in the subspace of independent densities. In 
particular, there is only one single elastic constant, which is A = 9 for the one- dimensional 
model. 

Bounds for the maximal entropy can be obtained as follows. The configurations with 
highest frequency of monominoes of type a are arrangements where all a's fill one sub- 
lattice. As the other vertices can independently be filled by monominoes of type b or c, 
we get a lower bound for the entropy: s max > log a/2. The one dimensional problem, at 
its point of maximal entropy, is a genuine Bernoulli type ensemble of entropy log 2, so we 
found 

0.3466 « log y/2 < s max < log 2 w 0.6931. (59) 

The given counting problem has first been solved by Baxter using Bethe's Ansatz. 
Introducing chemical potentials pi,fi2, A*3 fc> r the different monominoes, the grand-canonical 
potential reads M] 

/ / x , 8(z 1 z 2 z 3 )^ 

<PWl, f*2, f*3) = Wg 7= 7= W , 60 
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Figure 9: The entropy maximum {p = 0.50) 



where u,v,w are the positive real solutions of 

x(3B - xf = 4, u>v>w, (61) 

where B is given by 

Z\Z 2 + Z 2 Z 3 + Z1Z3 

3(^i z 2 z 3 ) 3 

It has to be stressed again that this solution was obtained by making use of periodic 
boundary conditions. The result agrees, on the other hand, with the corresponding free 



limit, as can be shown |56| by generalizing an argument of G. Kuperberg |4 



The value at the maximum entropy is 



~ log i w 0.4315, (63) 



which was already obtained before as the residual entropy of square ice [45 



It is instructive to analyse the entropy function along a line of reflection symmetry, 



see figure [L2|. Note that the entropy at its maximum is flat, indicating a vanishing 
second derivative, hence a phase transition at s — s max . The phase transition is of fluid- 
solid type: below the phase transition, monominoes arrange on the lattice homogeneously, 
whereas above the phase transition one sublattice is preferred in occupation ||52|| . The 
phase transition causes a square-root singularity in the covariance matrix M. 



In fact, explicit expressions for the entropy can be given in the two phases |p6[ . It 
turns out that both functions are analytic on the whole interval < p < | and intersect 
smoothly at pi — |. They can therefore be polynomially approximated at the entropy 
maximum. The absolute values of the coefficients coincide up to fourth order; with 
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Figure 10: Below the phase transition (p = 0.77) 



r = p — g, we find 

3 4 1 Q 2 1 3 

s ( r ) = _ log - - ^ — Irl 3 - 4r— r 4 + 0(r 5 ). (64) 
2 3 3! 2 4! 2 v ; v ; 

We conclude that the square lattice three-colouring model violates the second random 
tiling hypothesis! This means that the usual reasoning about elastic theory |3(| is not 
applicable here. 

It would be interesting to compare this result with an expansion of the entropy in 
terms of variables taken from the height representation of this model ||16|| . 

6.4 Hard hexagons — a monomino tiling with unexpected entropy 

Our last example belongs to a tiling class with two kinds of monominoes which we call 
particles and holes together with nearest neighbour exclusion for the particles. Its one- 
dimensional version has already been discussed above. The two-dimensional generaliza- 
tions to the square lattice and to the hexagonal lattice correspond to the ensembles of 
hard squares [ITDl |2T| and hard hexagons || [TTJ. Both models show a kind of fluid-solid 



phase transition: the low density phase is homogeneous, each sublattice has the same 
particle density, whereas at high density one sublattice is preferred in occupation. There 
is no obvious symmetry in these models, and therefore no preferred point in the phase 
space where the entropy maximum should be located. 

It is interesting to mention that the three-colouring model discussed above was in- 
terpreted as an 'approximating' hard squares model, with two types of holes [§, [52"1| . As 



the hard squares model itself has escaped exact solution so far, we concentrate on the 
hard hexagon model in the following. The solution of this model was given by Baxter 
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Figure 11: Above the phase transition (p = 0.90) 



|9|, |TlJ. This result is again obtained from periodic boundary conditions. It is easy to see, 
however, that it coincides with the result for the corresponding free boundary conditions, 
as each free rectangular patch can be transformed into a periodic one by adding a row of 
holes. 

A phase transition occurs at p c = ~ 0.2764 (the maximal density being normalized 

to |) with a critical exponent a = | in the covariance matrix. The entropy is a concave 
function with a quadratic maximum 

W ~ 0.333243 at p w 0.162433. (65) 

Note that the maximum occurs nearly at the point where the hexagons occupy half of the 
lattice area. The small deviation of the entropy from | lead to an incorrect prediction 
in 1978. The elastic constant turns out to be A ~ 24.0807. 

Joyce |36[ used the theory of modular functions to write the particle activity z as an 
algebraic function of the mean density p. Applying Joyce's results, one can in fact find 
closed expressions for s max ,p and A. As these are quite complicated, we omit them here 

m. 



7 Concluding remarks 

We have shown, both conceptually and by a number of examples, that the densities of tiles 
in random tiling ensembles provide natural parameters to investigate the phase diagram 
of such models. This is then independent of the existence of a height representation, but 
fully consistent with it if it happens to exist. 

10 For a short history of the model, see chapter 14.1 of Baxter's book [jLjJ . 
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Pi 

Figure 12: Entropy of the three-colouring model along a line of reflection symmetry (p2 = /O3) as 
a function of p\ (heavy curve) . The dotted curve represents the system with exclusion restricted 
to one dimension. 

Whereas the re-interpretation of exactly solvable models directly leads to a number of 
random tilings with analytic expressions for the entropy function, it is much more difficult 
to extract information about diffraction properties via Fourier transforms because this also 
requires the knowledge of the autocorrelation. For simple models such as the domino tiling 
and the rhombus tiling, it is possible || to compute the diffraction Bragg part explicitly 
and also the behaviour of the diffuse background in the thermodynamic limit. 

The next step in the analysis of random tilings is certainly the application of the con- 
cept developed so far to quasicrystalline random tilings. Although there are recent exact 
results on rectangle-triangle tilings |22J, [571, the interpretation of some of the results, for 
example the shape of the entropy curve in the solvable regime, is subject to controversial 
discussion. It may turn out that the framework of statistical mechanics proposed here 
gains additional insight into the phase space of the tilings. 
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Appendix: Height representation of the rhombus tiling 



Here, we briefly describe the connection between our description of the rhombus tiling to 
the description by Henley [^, |3(|. In particular, we show the connection between the 



elastic constants defined by both approaches. 



As was first pointed out in f\A \, each configuration of the rhombus tiling can be 
obtained by projecting a roof-type surface in the cubic lattice to the diagonal hyperplane. 
In this way, the rhombus tiling possesses a natural height representation. As the average 
slope of this surface determines the average densities of the different tiles and vice versa, 
we can parametrize the ensemble alternatively by some slope variables. 

Let us describe this more formally. Our vector space is R 3 with the standard scalar 
product < •, • >. The diagonal hyperplane in R 3 can be described by its normal vector 
e L = ^=(1,1,1)*. We call the space spanned by e 1 - perp(endicular) space V 1 - and the 

complementary hyperplane physical space Vk Let us denote the projections of vectors 
x E R 3 to these subspaces by x 1 - resp. x". The perp strain tensor E is a linear map 

E : V 11 -> V ± . (66) 

E defines a hyperplane in R 3 . By definition, E = corresponds to V". We define matrix 
elements Ei, E 2 of E: 

E x e{ = Ee\, = Ee l (67) 

A normal vector to this hyperplane is given by (1 — Ex, 1 — E 2 , 1 + E\ + E 2 y. 

We can easily relate the average slope of a surface parametrized by E = (E\, E 2 ) to 
the mean densities of the projected tiles: 

Pi = \{l-E 1 ), p 2 = l -{l-E 2 ). (68) 

The quadratic form ( [47|) can (with r = p — |) be rewritten in the form 

s 2 (E, E) = -±~.^.(E 2 1 + E X E 2 + E\). (69) 

This is the expression given by Henley [^, [3D| up to a constant of -7= which arises from 
different normalizations of unit vectors in physical and perp space in his setup. 
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